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We discuss collective excitations of a trapped dilute Fermi 
gas within a hydrodynamic approximation. Analytical results 
are derived for both high- and low-temperature limits and are 
applied to 40 K and 6 Li systems of current experimental inter- 
est. We identify spectral signatures which can be used to 
detect the onset of Fermi degeneracy. Also, we find an inter- 
esting class of internal excitations with an unusual spectrum. 
Some of our results are relevant to the case of trapped bosons 
as well. Our analysis suggests several experiments which ad- 
dress fundamental problems of collective motion in quantum 
fluids. 

A new era in low-temperature physics was heralded 
by the production in 1995 of Bose-Einstein condensation 
(BEC) in the alkali gases 87 Rb, 23 Na, and 7 Li. These 
dilute-gas systems are governed by atomic interactions 
that are very well characterized, so they provide systems 
in which quantum many-body physics can be explored 
from first principles. They have also been used for a re- 
markable variety of experiments on coherent matter wave 
generation and propagation, and have stimulated much 
experimental and theoretical work. 

Recently, impressive experimental results concern- 
ing the quantum degenerate regime of a dilute gas of 
trapped fermionic 40 K atoms have been presented |J 
thereby enabling fundamental investigations of quantum- 
degenerate Fermi systems with controlled particle inter- 
actions. Furthermore, reports of magneto-optical trap- 
ping of 6 Li have been presented recently. j| Several the- 
oretical results concerning the equilibrium properties of 
such a gas have been presented. [||-^) 

In the present paper, we examine the spectroscopy of 
collective modes of excitation of the Fermi gas. The study 
of this spectrum has been fruitful in studies of BEC. 
These modes correspond to sound waves in a homoge- 
neous gas, but for a confined gas their spectrum is dis- 
crete. They can be observed as shape oscillations of the 
atomic cloud, induced by perturbations of the trapping 
potential. Such measurements have been carried out for 
Bose gases, for cases of a nearly pure condensate |Io| , |lT| , 
for a partially-condensed gas at temperatures compara- 
ble to the BEC transition temperature T , fl^,|l|| and in 
the normal state |l3|]. For T < 0.6Xb the observed col- 
lective excitation frequencies are in good agreement with 
predictions of first-principles mean field theory |p^ , |l5| , 
and, for gases with more than a few thousand atoms, 
they agree well with the predictions of a simpler hydro- 
dynamic theory, |l6|,[l7j which resembles the approach we 
develop in this paper. 

As is well known jl8| , the collective modes of a quan- 
tum fluid admit simple descriptions in the collisionless 



and the hydrodynamic regimes. When the lifetime r of 
the quasiparticles is much longer than the characteristic 
period of motion (i.e. lot 3> 1 for atoms in a trap of fre- 
quency uj), there are few scattering events per sound os- 
cillation, and the restoring force is due the self-consistent 
mean field of the gas. Wave motion encountered in this 
limit is designated "zero sound". For the hydrodynamic 
regime lot <C 1, on the other hand, collisions ensure lo- 
cal thermodynamic equilibrium. To attain the hydro- 
dynamic regime in an ultracold gas of fermionic atoms, 
it is necessary to trap at least two different hyperfine 
states, since the (p-wave) interaction between atoms in 
the same hyperfine state is completely suppressed below 
T ~ 100/iK. Experimentally, trapping two hyper- 
fine states has been found to provide the mechanism for 
collisional rethermalization needed for evaporative cool- 
ing. f| This paper treats the case of a two-component 
Fermi gas in the hydrodynamic regime. We present a 
systematic analysis of the modes in the Fermi degenerate 
regime T < T F and in the classical limit T >0(T F ). Our 
analytical results for the density dependence of collective 
excitation frequencies suggest a means of detecting the 
onset of Fermi degeneracy in the gas at low T. We relate 
our results to current experiments on 40 K and 6 Li. 

The problem of collective modes of spatially confined 
two-component Fermi systems is well-known in the field 
of nuclear physics. However, due to the relatively few nu- 
cleons in a typical nucleus, these systems tend to be most 
appropriately described in the zero sound limit contrary 
to the situation described by the present paper. j^TJ 

In the limit of low energy/long wavelength excitations, 
one can treat the dynamics of a quantum gas by a semi- 
classical approximation, via the Boltzmann equation for 
the function /(r,p, t), which describes the distribution 
of particles of mass m with position r and momentum p. 
In the hydrodynamic regime, collisions ensure local equi- 
librium with a given mean velocity u(r, t), temperature 
T(r,t), and mass density p(r,t), since they conserve net 
momentum, energy and particle number. Assuming that 
the gas is otherwise ideal, we look for solutions to the 
Boltzmann equation that yield the Fermi-Dirac distribu- 
tion /(r,p,i) = {exp[/?(r,t)£(r,i)] + l}" 1 with £(r,t) = 
[p-mu(r,f)f/2m-fi F (r,i) and /3(r,t) = l/k B T(r,t). 
In this limit, the linearized momentum and energy con- 
servation laws are: 

p (r)9 t u(r, t) = -VP(r, t) + p(r, p, i)f (1) 



^V[P (r)u(r,t)] + ^ (r)u(r,<)f. (2) 



d t P{r,t) : 

with f = — VF poi (r)/m. Here, V pot (r) is the trap- 
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ping potential and p(r,t) = po(r) + 8p(r,t), P(r,t) = 
Po(r) + SP(r,t) with po( r ) an d Po(r) being the equi- 
librium mass density and pressure respectively. As we 
are interested in modes for which the densities of the 
two states oscillate in phase, we have defined the to- 
tal density p(r,t) = £ ff f d 3 k / (2n) 3 f (r , p,i) = £ ff(0ff (r) 
where denotes the sum over the two trapped hyper- 
fine states. Likewise, P(r, t) and u(r, t) denote the total 
pressure and mean velocity associated with both species. 
We take the number of atoms in each hyperfine state to 
be the same, so that there is a single chemical poten- 
tial and p(r) = 2p a (r). The "spin" -modes, for which the 
two components of the gas do not move in phase, will 
be purely diffusive in the hydrodynamic regime, as the 
collisions suppress any difference in the local velocity of 
the two components. |21| Using the continuity equation, 
d t p(r,t) = -V[p(r,<)u(r,t)], Eqs. (@)-(§) can be written 
as a closed equation for the velocity field pa] : 



5Po 
3po 



V(Vu)+V(uf) + -f(Vu). 



(3) 



Hereafter, we take the trap to be an isotropic harmonic 
oscillator, i.e. V po t = muj 2 r 2 /2. The hydrodynamic equa- 
tions are then solved in several limiting cases. We first 
consider the limit T -C Tf , where the pressure is a func- 
tion only of the density and we invoke the semiclassical 
(Thomas-Fermi) approximation to the equilibrium den- 
sity distribution, 



Po(r) 



8Nm 



I 2 2\ 3 / 2 



(4) 



where N is the total number of atoms in the trap, and 
r F = ^2/ip/mw 2 , with p F = (3 AT) 1/3 From Eq.(§) 
we obtain Po/po = uj 2 (rp — r 2 )/5, so Eq.(0) becomes 



dfu = — V [(r| - r 2 )Vu] - w 2 V(ur). (5) 

o 

As the right hand side of this equation is a pure gradient, 
the sound modes can be found by considering potential 
flow. Equation (|lf) and the continuity equation then de- 
couple, to yield 



d 2 8p 



5* + |fc 



8p, 



(6) 



Eq.(^|) clearly separates in spherical polar coordinates 
(r, 9, </>), so we invoke the eigenfunction expansion, 



singular points, and the condition that 8p„i(r) be finite 
at all r yields the eigenfrequencies |23] 



2uj 



uj s = —= \/n 2 + 2n + In + 3Z/4 
v3 

and the (unnormalized) eigenfunctions 



(8) 



Spni (r) = r l y/rl-r 2 2 F X [-n, I + 2 + n, / + 3/2; (r/r F ) 2 } 

(9) 

where n = 0,1,2,..., and 2 Pi is the usual hypergeometric 
function [B4J, which is equivalent to the Jacobi polyno- 
mial p^ 1 / 2 ' 1 / 2 ) , For I = 0, the solutions reduce to 

8p(r)n0 oc sin[(2n + 2) arccos(r/rF)]/7" 

cx yJl-r 2 /r 2 F U 2n+1 (r/r F )/r (10) 



with U n (x) being the Chebyshev polynomials. Note that 
the zero-frequency solution, ujqq = 0, is associated with 
a mode, 8p 00 (r) = (1 — r 2 /^) 1 / 2 , which corresponds 
to an infinitesimal change in the number of particles in 
the system [see Eq.(Q)]. It is therefore not a physical 
sound mode. All other modes described by Eqs.(|3||9|) are 
particle conserving [i.e. J drr8p n (r) — 0], as required. 
The lowest frequency for a given I is ujiq — \ftuj, a re- 
sult first derived by Griffin et al. p3| The u>\q = uj 
mode corresponds to the lowest center-of-mass oscilla- 
tion. As noted by Griffin et al. [Q , only modes for which 
Vu ^ constant depend on the statistics (i.e. Po/po) of 
the atoms as can be seen directly from Eq.(||). By cal- 
culating the complete sound spectrum, we have found 
these modes. From V 2 [r'Yj m (#, <j))} = 0, we see that for 
any angular momentum, the lowest observed mode does 
not depend on the statistics of the gas whereas the higher 
ones (n > 1 for I = and n > for I > 0) do. 

We now consider the high T regime where the de 
Broglie wavelength is much smaller than the interpar- 
ticle spacing and the gas approaches the classical limit. 
However, we assume sufficiently low temperatures such 
that the density is high enough for the hydrodynamic 
approximation to be valid. To find the sound modes in 
this regime, we need to return to Eq.(^|) as the pres- 
sure is now not related to the density in a simple way. 
For I = 0, the velocity field is irrotational and we write 
u = V[p n (x)], where p n is a polynomial in x = r/r^p 
with ttf = \J k-gT /miu 2 . Using Po/po = k^T/m, we 
find the sound frequencies 



Sp(r,t) = Y / e iUnlt Spni(r)Y lm (9, 

nlm 



(7) 



where n is an eigenvalue index, u n i a sound eigenfre- 
quency, and 8p n i (r) the corresponding density fluctuation 
eigenfunction (with a normalization to be determined by 
boundary conditions). Substituting Eq.(0) into Eq.(g) 
gives an ordinary differential equation with three regular 



uj 0n — uy/lOn/3 + 4 



(11) 



with n = 0, 1, 2, . . .. The corresponding density fluctua- 
tions are given by 

8p(r) oc (2/r + d r )[e'^ 2r ^ 2 d r pn(r/r F )} (12) 

where d x p n (x) = x[b + b 2 x 2 + . . . + b 2n x 2n ]- As ex- 
pected, we recover the u> s = 2u mode which, since 
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Vu = constant, is independent of the particular form 
of Pq/pq. However, the higher modes (Vu ^ constant) 
do depend on the statistics as can be seen by comparing 
Eq.(|[) with I = and Eq. (|TT|) . For general angular mo- 
mentum I, the analysis is somewhat complicated by the 
fact that the velocity modes are not irrotational. Writing 
u = -VV> + V x A and 2rVu/3 = -V^ + V x A, we 
obtain from Eq.(|^) the coupled equations: 



%4> 

d 2 (v x a) 



5k B T 



V 2 V>-w 2 (VV>- V x A)r -lu 2 iP (13) 



3m 
-w 2 V x A 



(14) 



They are solved by writing -0„; = p n i(x)x l Yi m (9, c/>) with 
Vni = clqi + a,2ix 2 + . . . + a2 n ix 2n . The spectrum is 



-(Z + 2n+2/5) 



± yj(l + 2n + 2/5) 2 25/9 - 1(1 + 1)8/3 



(15) 



with n = 0,l,2,.... The result u> s = y/luj for the lowest 
mode for a given I independent of statistics is recovered 
as expected. Again, by comparing Eq.(|) and Eq.© 
we see that the higher modes depend on the functional 
form of the density profile. The corresponding density 
fluctuations obtained from the continuity equation are 
Sp n i oc exp(-x 2 /2)x l [a i + a 2 ix 2 + . . . + a 2n ix 2n }Yi m (9 , (p) . 
Using this, it is easy to show that the density fluctuations 
are particle conserving. The high T sound frequencies 
are independent of T which is a special property of the 
harmonic trap. One should note that the hydrodynamic 
approximation starts to break down for the high lying 
modes (n large) when the characteristic wavelength of 
these modes becomes comparable to the mean free path. 

The eigenfrequencies given by subtracting the square 
root in Eq.([l5]) (for n > 1) represent an interesting class 
of solutions. These modes have decreasing frequency with 
increasing wave number which is characteristic for the 
so-called internal waves. Internal waves are present in 
general when the equilibrium density profile is stratified 
by an external force. They are, contrary to ordinary 
sound, essentially driven by the external force giving rise 
to the unusual relation between frequency and wave num- 
ber. 25 A detection of these modes would be a very 
useful experimental verification of a general result in the 
theory of hydrodynamics. 

We will now turn to possible experimental implica- 
tions of the results presented. As mentioned above, in 
order to achieve the hydrodynamic limit lot <C 1 one 
should trap the relevant fermionic atoms in two hyper- 
fine states. We need to examine the effect of the Fermi 
blocking on the scattering rate of the atoms. In the 
presence of a Fermi sea, the low energy lifetime r of a 
quasiparticle for low temperatures is T c i ass /T oc (T /Tp) 2 
where t~^ ss ~ p a ao(v)/m is the classical scattering rate, 
Tp = pp/kB, and (v) is the thermal velocity of the 



atoms. For the energies relevant for the present prob- 
lem, only s— wave scattering between atoms in different 
hyperfine states is important and o~q = Aira 2 where a 
is the s— wave scattering length. The (T/Tp) 2 factor 
reflects the increasing Fermi blocking of the scattering 
states with decreasing temperature. |L^] For T > 0(Tp) 7 
the Fermi blocking becomes negligible, and r _1 ~ T cilss- 
As (v) cx T 1 / 2 and p(0) cx T" 3 / 2 , we have t^ ss oc T" 1 . 
Thus, both for T — ► and T — > oo the gas is in the colli- 
sionless regime and we need to examine in which region 
of the phase diagram the hydrodynamic description pre- 
sented in this paper is valid. We will concentrate on the 
two atoms 6 Li and 40 K and assume a trapping frequency 
of v = lo/2it = 144Hz approximating typical experimen- 
tal conditions. Of course, by altering the trapping fre- 
quency one can change the regions where the hydrody- 
namic approximation is valid. 

For 6 Li, the scattering length is a ~ — 2160ao where ao 
is the Bohr radius. |2(J This large value should mean that 
the hydrodynamic description is valid in large regions of 
the phase diagram. Indeed, using v — 144Hz and assum- 
ing that there is 10 6 particles trapped in each hyperfine 
state, we obtain that the hydrodynamic approximation is 
accurate down to temperatures T ~ O.lTp. For such low 
temperatures, the gas is in the Fermi degenerate regime 
and we expect Eqs.(|^-[)|) to give a good description of 
the sound modes. As T is increased to T ~ Tp, Po/po 
approaches the classical form and the sound spectrum is 
given by Eq. (O) . Since the scattering length for 6 Li is 
so large, the collisionless limit is only reached for a rel- 
atively low density ( T > Tp for N=2 x 10 6 ). Thus, 
for a wide temperature range the hydrodynamic approx- 
imation is valid and the sound spectrum should be inde- 
pendent of T and given by Eq.([15|). The onset of Fermi 
degeneracy for T < Tp should show up in a T depen- 
dence of the sound spectrum as it changes from Eq.(|T5|) 
to Eq.(^) with decreasing T reflecting the change in the 
functional form of Pq/pq. Eventually, for T 3> Tp we 
have lot ~ 1 and the damping of the modes becomes 
large ^8[?,^7). For higher T, the gas is in the collision- 
less regime. In this regime, the frequencies approach the 
single particle spectrum, e.g. lo s ~ 2nu for I = 0, as the 
effect of the interactions decrease with decreasing den- 
sity pl[] . Of course, by trapping more/fewer atoms one 
can increase/decrease the temperature region where the 
hydrodynamic description is valid. 

The spectrum given by Eq.(^) coincide with the T = 
collective mode spectrum in the BCS state. |^8| Thus, the 
density-fluctuation spectrum in the hydrodynamic limit 
of the normal phase and in the superfluid phase is the 
same. This generalizes the well-known result for homo- 
geneous systems, to the case of trapped systems. The 
effect of the superfluid correlations on the various collec- 
tive modes in general has been examined in detail for 
homogeneous systems. 
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For 40 K, the situation is very different since the scat- 
tering length 152ao < a < 542ao is much smaller. M] 
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Therefore, in order to reach the hydrodynamic regime 
one needs to trap at least ~ 10 8 atoms in each hyperfine 
state. Alternatively, one could increase the scattering 
rate by adjusting an external magnetic field to achieve 
a Feshbach resonance. This latter procedure is also ex- 
pected to be necessary to obtain a large and negative 
scattering length in order to get an experimentally re- 
alistic transition temperature T c for the BCS transition. 
Therefore, although the treatment of the sound modes 
put forward in this paper does not describe the present 
experimental situation for 40 K, one might expect it to be 
relevant for future experiments. 

Of course, Eq.(fLll) and Eq.(15) are still valid for high 



temperatures for a bosonic gas in the hydrodynamic 
regime. As an additional result, we therefore suggest 
that one should also be able to observe the transition 
from the classical to the quantum regime for a Bose gas 
by looking at the change in the sound spectrum. The 
modes for high T are presented in this paper whereas the 
modes in the BEC limit are described in ref. jl(| . 

For an anisotropic gas, a calculation of the higher 
(Vu ^ constant) modes is somewhat cumbersome al- 
gebraically. However, although the spectrum will be dif- 
ferent, the qualitative behavior is the same as for the 
isotropic case presented here: the modes are still inde- 
pendent of T in the classical limit and one should be able 
to detect the onset of Fermi degeneracy by looking for a 
T-dependence of the modes. Thus, the main conclusions 
of the paper remain valid for the anisotropic case. 

In conclusion, we have analyzed the collective modes 
of a dilute trapped Fermi gas in the hydrodynamic limit. 
Analytical results both in the low and high T limit for 
modes which depend on the density profile of the cloud 
were presented. This lead to the proposal of a straight- 
forward method to detect the onset of Fermi degeneracy. 
We found an interesting class of internal modes which 
have decreasing frequency for increasing wave number. 
Furthermore, the analysis presented for the high T spec- 
trum should also be valid for the case of trapped bosons. 
Our analysis should have direct relevance to current ex- 
periments on 6 Li atoms and to possible future experi- 
ments on 40 K. 

We appreciate useful discussions with D. J. Allwright, 
H. Ockendon, and D. S. Jin. 
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